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Yields, correlation shapes, and mean transverse momenta pr of charged particles associated with 
intermediate to high-pT trigger particles (2.5 < pr < 10 GcV/c) in d+Au and Au+Au collisions 
a t v^nn" = 200 GeV are presented. For associated particles at higher pr > 2.5 GeV/c, narrow 
correlation peaks are seen in d+Au and Au+Au, indicating that the main production mechanism 
is jet fragmentation. At lower associated particle pr < 2 GeV/c, a large enhancement of the near- 
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(A<f> ~ 0) and away-side (A(/> ~ 7r) associated yields is found, together with a strong broadening 
of the away-side azimuthal distributions in Au+Au collisions compared to d+Au measurements, 
suggesting that other particle production mechanisms play a role. This is further supported by the 
observed significant softening of the away-side associated particle yield distribution at A(j> ~ ir in 
central Au+Au collisions. 

PACS numbers: 25.75.Nq, 25.75.Bh 



I. INTRODUCTION 

The goal of ultra-relativistic heavy-ion collisions is to 
create a system of deconfmed quarks and gluons at high 
temperature and density and study its properties. In the 
initial stage of the collision, hard scatterings between par- 
tons in the incoming nuclei produce high transverse mo- 
mentum (pt) partons that fragment into jets of hadrons 
with a clear back-to-back di-jet signature [l[ . In Au+Au 
collisions, hard partons traverse the hot and dense col- 
ored medium, thus probing the medium through energy 
loss 0-13. 

In-medium jet energy loss was first observed at the Rel- 
ativistic Heavy Ion Collider (RHIC) as a suppression of 
hadron spectra at high pr [1,0] in Au+Au collisions with 
respect to p+p collisions. The jet-like structure of hadron 
production at high px was confirmed by measurements 
of the azimuthal angle difference Acf> distributions of as- 
sociated particles in a certain range of px with respect to 
a trigger hadron at a higher px [l| . At the highest px , a 
suppression of the away-side yield (around Acj) ~ tt with 
respect to the trigger particle) by a factor 3-5 is observed 
0] . This suppression is consistent with theoretical calcu- 
lations that incorporate in-medium energy loss [1, Q . At 
lower of the associated particles, a strongly broadened 
away-side structure is seen in Au+Au collisions, and the 
associated yields on both the near-side (Acj) ~ 0) and 
away-side (A<j> ~ n) are enhanced [T3, EH ■ A number of 
possible explanations of the away-side broadening at in- 
termediate px have been put forward, ran ging from frag- 
mentation products of radiated gluons [l2T,ll3j| to medium 
resp onse and the possibility of a Mach-cone shock wave 

The p^ 19 range used in previous studies [ToL ITl| is 
the region where the p/-7r ratio is large. The large 
baryon/meson ratio has been interpreted as being due to 
coalescence/recombination of quarks, which could also 
have an impact on the jet-like correlation yields, espe- 
cially for trigger particles in the px range 2.0 to 4.0 
GeV/c where coalescence/recombination products (l8l - 
|20| may be present. 

In this paper, we present a systematic exploration of 
the azimuthal di-hadron correlation shapes and yields 
with centrality and px of the trigger (Px~ 19 ) and asso- 
ciated hadrons (p T ssoc ), to investigate the change from 
broadened correlation peaks with enhanced yields at low 



Px to suppressed away-side yields at high px- The analy- 
sis is performed on the large statistics sample of Au+Au 
collisions at -y/sjvw = 200 GcV collected by the STAR 
experiment in the RHIC run in 2004. The d+Au data 
sample from the year 2003 is used as a reference where no 
hot and dense matter is formed, because the minimum 
bias p+p data collected by STAR has limited statistics. 
Earlier studies have shown that di-hadron correlations in 
p+p and d+Au collisions are similar 



II. EXPERIMENTAL SETUP AND DATA SETS 

The measurements presented in this paper were per- 
formed with the STAR detector at RHIC [2lj]. Charged 
tracks are reconstructed with the Time Projection Cham- 
ber (tpc) [13. 

For Au+Au collisions, two different online event selec- 
tions (triggers), minimum bias and central, were used. 
The central trigger selection was based on the energy 
deposited in the two Zero-Degree Calorimeters (ZDCs) 
which measure spectator fragments and small-angle par- 
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ticle production |21[. The trigger selected the most cen- 
tral 12% of the total hadronic cross section, based on a 
maximum energy deposited in the ZDCs and a minimum 
multiplicity in the Central Trigger Barrel (CTB) [23|. 
The central trigger also uses time information from the 
Beam-Beam Counters (BBC) to restrict the primary ver- 
tex position z vtx to be within approximately ±30 cm of 
the center of the detector along the beam direction. The 
minimum bias (MB) trigger is based on a ZDC coinci- 
dence (a threshold amount of energy in each ZDC) and 
required a minimum multiplicity in the CTB to reject 
non-hadronic interactions. For the minimum bias sam- 
ple, events were selected to have \z vtx \ < 25 cm. A total 
of 21M minimum bias events and 18M central triggered 
Au+Au events were used. 

For d+Au collisions, the minimum bias trigger was de- 
fined by requiring that at least one beam-rapidity neu- 
tron impinge on the ZDC in the Au beam direction. 
The measured minimum bias cross section amounts to 
95 ± 3% of the total d+Au geometric cross section pj. 
For d+Au events, the distribution of primary vertices 
along the beamline was wider than during the Au+Au 
run. The events were selected to be within ±50 cm from 
the center of the detector along the beam line. A total 
of 3.4M d+Au events were selected for this analysis. 

The Au+Au events are further divided into centrality 
classes based on the uncorrected charged particle multi- 
plicity in the range \rj\ < 0.5 as measured by the TPC. We 



4 



present results for the following centrality ranges: 0-12% 
(from the central triggered data set), 20-40%, 40-60%, 
and 60-80% (from the MB data set) of the total hadronic 
cross section, with 0% referring to the most central col- 
lisions. 



III. DATA ANALYSIS 

Di-hadron correlations are constructed using charged 
particles measured in the TPC. All particles are selected 
to have pseudo-rapidity in the range —1.0 < r\ < 1.0, so 
that they fall well within the TPC acceptance. To reject 
background tracks at high p^, tracks were required to 
have at least 20 measured points in the TPC (out of 45) 
and a distance of closest approach (dca) to the event 
vertex of less than 1 cm to reduce the contribution from 
secondary particles. 

The results are corrected for single particle acceptance 
and detection efficiency as well as for the pair acceptance 
as a function of A(f>. The single particle reconstruction 
efficiency as a function of tj, px, and centrality is deter- 
mined using hits from a Monte Carlo simulation which 
are embedded into real data events. The tracking effi- 
ciency depends sensitively on the gain in the proportional 
readout chamber of the TPC and thus on the atmospheric 
pressure. Uncertainties in details of these effects give rise 
to an overall 5% systematic uncertainty in the absolute 
yields given in this paper. In most cases, the uncertainty 
from the background subtraction as described in the next 
section is larger than the systematic uncertainty from the 
tracking efficiency. The TPC sector boundaries introduce 
a dependence of the pair acceptance on angle difference 
Acf>, which was determined from mixed events. No Arj 
pair acceptance correction has been applied. A small in- 
efficiency due to tracks crossing inside the TPC volume 
affects the asssociated hadron distribution at small pair 
separation in (A77, A(j>). This effect manifests itself as a 
reduced efficiency for small but finite A0, at positive or 
negative A<p, depending on the sign of the curvature of 
the associate track. A correction was performed by first 
curvature-sorting the distributions and then reflecting a 
few bins from the unaffected area to the area where the 
inefficiency occurs, thus restoring the symmetry between 
positive and negative A<f>. 

Figure Q] shows an example azimuthal angle difference 
distribution for trigger particles with 4.0 < p^ lg < 6.0 
GeV/c and associated particles with 1.0 < p T ssoc < 2.5 
GeV/c in 0-12% central Au+Au collisions. The distribu- 
tion is divided by the number of trigger hadrons to give 
the associated yield per trigger hadron. The associated 
hadron distribution contains a background of uncorre- 
cted particles which has a cos (2A<j>) modulation due to 
the correlation of all particles with the reaction plane 
through elliptic flow, v 2 - We model the background us- 
ing the function B(l + 2(v t 2 ri9 )(v% ssoc ) cos (2A(f>)), where 
the v 2 values are from separate flow measurements [24j |. 
The function is normalized to the data in the region 
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FIG. 1: (Color online) Azimuthal distribution of associated 
charged particles with 1.0 < p^ saoc < 2.5 GeV/c with respect 
to trigger particles with 4.0 < p% ia < 6.0 GeV/c in 0-12% 
central Au+Au collisions. The curve shows the modulation 
of the background due to elliptic flow V2 and the grey band in- 
dicates the uncertainty on the elleptic flow of the background 
(see text). 



0.8 < \Acf)\ < 1.2, where the signal is apparently small, 
and then subtracted. This background normalisation 
procedure is often referred to as the ZYA1 ('Zero Yield at 
1 radian') or ZYAM ('Zero Yield At Minimum') method 

The method of normalising the combinatorial back- 
ground level in the region around \A(j)\ = 1 has first been 
used for di-hadron correlations at higher momenta (7l[l0j. 
where there are narrow peaks on the near- and away-side, 
separated by a largely 'signal-free' region. At lower pt, 
the correlation peaks are broader, and there is no clear 
signal-free region, so that the background normalisation 
is more ambiguous. In addition, due to the larger combi- 
natorial background, the elliptic flow modulation of the 
background is of similar size as the trigger-associated- 
hadron correlation signal. The ZYAl method provides 
a simple prescription to separate signal and background, 
which we will use throughout the paper. An alterna- 
tive approach would be to decompose the correlation 
shape using a fit function with components represent- 
ing the flow modulation of the background and the as- 
sumed shapes of near- and away-side correlation peaks 
[26l [27| . The unsubtracted azimuthal hadron distribu- 
tions are provided in the Appendix and can be used for 
such a procedure. 

The nominal v 2 value used for the subtraction is the 
mean of the v 2 measured using the reaction plane method 
with the Forward TPC and the four-particle cumulant 
method, which have different sensitivity to non-flow ef- 
fects and flow fluctuations (line in Fig. [1]). The dif- 
ference between the two results is used as the estimate of 
the systematic uncertainty in v 2 and this range is shown 
by the band in Fig. [T] For the d+Au results a constant 
pedestal (normalized in the same A<f> range) was sub- 
tracted. 
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IV. RESULTS 
A. Azimuthal di-hadron distributions 

Fig. [5] shows the background subtracted associated 
hadron A<j> distributions with 1.0 < p^ soc < 2.5 GeV/c 
for 4 centrality selections, 60-80%, 40-60%, 20-40% and 
0-12%, and 4 trigger selections, 2.5 < p% ig < 3.0 GcV/c, 
3.0 < V % ia < 4.0 GeV/c, 4.0 < p% 19 < 6.0 GcV/c, and 
6.0 < p^ 19 < 10.0 GeV/c. Results are presented for two 
different ranges in the pseudo-rapidity difference between 
the trigger and associated particles \Ar)\. The shapes are 
very similar for both Arj selections in all panels (there 
is an overall reduction in the away-side yields due to the 
smaller acceptance for |Ar/| < 0.7). For reference, the 
di-hadron distributions without background subtraction 
are shown in the Appendix, where also the V2 values and 
background normalization values (B) used to subtract 
the background are given. The systematic uncertainties 
on the V2 values for the background are shown by the 
bands around the data points. The d+Au results (open 
circles) are also shown for reference. 

In Fig. [2 top row, one observes that the jet-like correla- 
tions in peripheral (60-80% centrality) Au+Au collisions 
are very similar to the d+Au result, indicating that such 
correlations in peripheral Au+Au collisions can be de- 
scribed as a superposition of independent p+p collisions. 
The near- and away-side yields of associated particles in- 
crease with Pt' 9 ', as expected from parton fragmentation. 

For more central events, a significant increase of both 
the near- and the away-side yields is seen in Au+Au col- 
lisions relative to d+Au. The relative increase of the 
near-side yield is larger for lower p^ 19 (top row) than for 
higher p" 3 . For peripheral events, the near-side results 
for |Ar/| < 0.7 do not differ significantly from the full ac- 
ceptance results, demonstrating that the correlated yield 
is at relatively small A77, as expected from jet fragmen- 
tation. For more central collisions, on the other hand, 
a significant fraction of the associated yield is at large 
I At/ I > 0.7 for the lower indicating a significant 

long-range correlation in A77, possibly due to an inter- 
play between the soft bulk dynamics of longitudinal flow 
and jet-like di-hadron structure (28|. It has also been ar- 
gued recently that this long-range correlation in A 77 could 
be caused by long-range structures in the medium, due 
to density fluctuations in the medium (29j or color flux 
tubes [30l - l32l ] . The 'ridge'-like correlation structure in Arj 
is further explored in other STAR publications (2(| HH? 

It is interesting to note that the largest relative en- 
hancement of the near-side yield is observed for the lower 
p% 19 , 2.5 - 4.0 GeV/c. It has been suggested that parti- 
cle production in this momentum range has a large con- 
tribution from coalescence of quarks from bulk partonic 
matter [l8l - [20| . This production mechanism would not 
lead to jet-like structures. Trigger hadrons formed by this 
mechanism would increase the number of trigger hadrons, 



without increasing the associated yield, leading to a re- 
duced per-trigger associated yield, in contrast to what 
is observed in Fig. [2] The increased associated yield 
at intermediate pt indicates that if coalescence is a sig- 
nificant source of hadron production at intermediate pt, 
it has to generate angular correlation structure, either 
through shower-thermal coalescence [1^] or local fluctua- 
tions in the medium density or temperature, e.g. due to 
heating of the medium by the passage of a parton [36j . 
So far, most calculations of such effects are qualitative at 
best. Quantitative predictions for these processes should 
be made and compared to the data. Measurements with 
identified baryons and mesons as trigger and associated 
particles [37], HH can be used to explore the possible con- 
tributions from coalescence. 

The away-side yield of associated particles at low p 1 ^ 19 
(top row Fig. [2]) evolves significantly in both shape and 
yield with centrality: the shape becomes much broader 
than the d+Au reference and the yield increases. For 
20-40% central collisions the distribution becomes flat 
or slightly double-peaked, with a shallow minimum at 
A<j) = 7r. In the most central collisions the distribution is 
double-peaked for the lowest p^ 19 . With increasing p% 19 1 
the away-side shape becomes flatter. Overall, there is a 
smooth evolution of the peak shape with centrality and 
ptrig ^ rpj^g va j ue f f or w hich the away-side becomes 
flat or double-peaked decreases with centrality. Note that 
the double-peak shape is not seen in the raw signal in 
Fig. [U and only appears after the subtraction of the V2- 
modulated background. In that sense, the double-peak 
structure is generated by imposing a separation between 
flow and non-flow in the analysis of azimuthal correla- 
tions. This separation is not unambiguous and remains 
under active investigation. 

For the most central collisions, the broadening of the 
away-side structure is so large that the near- and away- 
side peaks may overlap, making it impossible to unam- 
biguously distinguish the correlation structure from the 
background without other inputs. For the present analy- 
sis, we have chosen to use the same background normal- 
ization procedure for all centrality bins and px bins, i.e. 
to normalize the i^-modulated background to the signal 
in the range 0.8 < \A(j>\ < 1.2 and subtract it. 

In Fig. [3] we focus on central data where the largest 
modifications of the correlation shapes and yields are 
found. The figure shows the correlation shapes in the 0- 
12% central event sample for different selections of p^? soc 
and p^ 19 . As in Fig. [21 results are given for the full Arj- 
acceptance (solid circles) as well as a restricted range 
|A?/| < 0.7 (squares), and for d+Au collisions (open 
circles). The distributions before background subtrac- 
tions and the background normalization and V2 values 
are given in the Appendix. 

On the near-side, we observe again a large increase 
of the yield in central Au+Au collisions compared to 
d+Au collisions. The yield depends on the Ar/-selection 
used, indicating that there is significant associated yield 
at Ar] > 0.7. The relative size of the enhancement de- 
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FIG. 2: (Color online) Background-subtracted azimuthal angle difference distributions for associated particles with pr between 
1.0 and 2.5 GeV/c and for different ranges of trigger particle pr, ranging from 2.5 — 3.0 GeV/c (left column) to 6 — 10 GeV/c 
(right column). Results are shown for Au+Au collisions (solid circles) with different centrality (rows) and d+Au reference 
results (open circles). The rapidity range is \r]\ < 1 and as a result the rapidity-difference |Ar;| < 2. Open red squares show 
results for a restricted acceptance of \Ar)\ < 0.7, using tracks within \r)\ < 1. The solid and dashed histograms show the upper 
and lower range of the systematic uncertainty due to the v% modulation of the subtracted background. 



pends on p9j? soc and p^ lg . The measured jet-like yield 
in d+Au collisions increases faster with p!^ 19 (going from 
left to right in Fig. [3]) than in Au+Au collisions, reduc- 
ing the relative size of the enhancement in Au+Au. The 
associated yield decreases with increasing pjf soc for both 
d+Au and Au+Au collisions, but the decrease is stronger 
in Au+Au, so that the measured yields in Au+Au ap- 



proach the d+Au results at the highest p^ ssoc . A sum- 
mary of the yields is presented in Fig. [5] (Section II V Cp . 

On the away-side, we observe a broadening and en- 
hancement of the yield in Au+Au compared to d+Au, 
except at 2.5 < p°h ssoc < 4.0 GeV/c (bottom row of 
Fig. [3]), where a broadening is seen, while the yield is 
smaller than in d+Au. For 6 < p l rf lg < 10 GeV/c (right- 
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FIG. 3: (Color online) Background-subtracted azimuthal angle difference distributions for different p^T 19 (columns) and p^ ssoc 
(rows) in 0-12% central Au+Au collisions (solid circles) and d+Au reference results (open circles). The rapidity range is \rj\ < 1 
and as a result the rapidity-difference |Ar;| < 2. Open red squares show results for a restricted acceptance of |A?7| < 0.7, using 
tracks with \rj\ < 1. The solid and dashed histograms show the upper and lower range of the systematic uncertainty due to the 
V2 modulation of the subtracted background. 



most column in Fig. [3]), a narrow peak appears at large 
passoc - m Au+Au, similar to what is seen in d+Au colli- 
sions and at higher px in Au+Au collisions Q ■ 

Although the shape of the away-side distribution 
changes with p^ 19 and p T ssoc , there seems to be no grad- 
ual broadening as a function of px'. the rising flanks of 
the away-side distribution are at similar A0 in the entire 
range 0.5 < p T ssoc < 2.5 GcV/c and 2.5 < p% 19 < 6. In 



fact, it could be argued that the away-side distribution 
is as broad as possible; there is no Atj> region without 
correlation signal. 

The broad away-side correlation structure in Au+Au 
collisions is a truly remarkable observation. Although 
some broadening of the away-side correlation in Au+Au 
collisions would be expected due to increased acopla- 
narity (kx) due to multiple scattering of the parton in 
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FIG. 4: (Color online) Background-subtracted associated hadron distribution in 0-12% central Au+Au collisions with 0.8 < 
p assoc < L0 GeV y c and 4 < ptrig < 6 GeV / c The co i or ed 

curves indicate fits with 3 different functional forms for the 
away-side shape: (left panel) symmetric identical Gaussian distributions (2 Gauss), (middle panel) adding a cosA(f> background 
distribution (2 Gauss + Cosine), and (right panel) adding a third Gaussian as the away-side jet-like component (2+1 Gauss). 
The histograms indicate the uncertainty on the background shape due to elliptic flow. 



the medium, the structures seen in Fig. [3] are broader 
than would be expected from such a mechanism [ljj . It 
has, however, been pointed out that kinematic selection 
effects on in-medium gluon radiation may lead to a non- 
trivial structure in the angular distributions [13|. It has 
also been argued that a fast parton may generate sound 
waves in the bulk quark-gluon matter which would lead to 
a Mach-cone shock wave Evidence for a conical 

emission pattern has been found in three-particle corre- 
lation measurements (39[. The broad structure seen in 
the di-hadron distribution could then be the projection 
of the conical pattern on A(f>. Another mechanism that 
may produce conical emission from a fast parton is QCD 
Cherenkov radiation El| . There are two other cal- 
culations that show a broad, double-peaked away-side 
structure without implementing a specific mechanism for 
conical emission: one is a 3D hydrodynamical calcula- 
tion which includes local density fluctuations in the ini- 
tial state [29[ and the other is the AMPT model [H[. It is 
worth noting that in the 3D hydrodynamical model, there 
is also no explicit introduction of hard partons or jets; 
the correlation arises purely from the medium. All these 
models should be confronted with the data presented in 
Figs and [3] as well as the three-particle correlation data 
in [39j. 

In addition to the change of the correlation shapes, a 
significant increase of the yields is seen in Au+Au col- 
lisions relative to d+Au collisions, for most of the px~ 
selections in Fig. [31 both on the near- and the away-side. 
The yield increase implies that trigger hadrons in Au+Au 
collisions are accompanied by a larger energy flow than 
trigger hadrons with the same transverse momentum in 
elementary collisions. This would be compatible with 
a scenario where the leading hadron is softened due to 
energy loss so that trigger hadrons in Au+Au collisions 
select a larger initial parton energy than in d+Au or p+p 
collisions. 



B. Away-side shapes 

To further characterize the broad shape of the away- 
side associated hadron distributions in Figs [2] and [3l the 
data were fitted with different parameterizations. Three 
different functional forms were used, all of which are 
based on the assumption that there is significant yield 
in a cone in (Arj,A(f)) around the away-side parton. The 
projection of this cone on A<p would give rise to two peaks 
symmetric around A(f> = tt. Fig. U shows the associ- 
ated hadron distribution after background subtraction 
for 4 < p% lg < 6 GcV/c and 0.8 < p^ ssoc < 1.0 GeV/c 
fitted with three different functional forms that include 
two Gaussian distributions at A(f> = tt ± A on the away- 
side. The left panel of Fig. 2] shows the simplest ansatz, 
using just two Gaussian peaks on the away side and a 
single Gaussian peak on the near side. To account for 
correlations induced by momentum conservation or rem- 
nant jet structure, we add a cos(A0) distribution (middle 
panel of Fig. [4| or a third Gaussian peak (with a different 
width and amplitude) at A</> = tt (right panel Fig. 
The three parameterizations are referred to as 2 Gauss, 
2 Gauss + Cosine and 2+1 Gauss, respectively. The best 
fit is obtained when using the 2 Gauss + Cosine function, 
as can be seen from the \ 2 values given in the figure. 

Figure [5] shows A, the angle (in radians) between the 
Gaussian peaks and A0 = tt, as a function of p°h ssoc , 
for the three different paramctrizations of the away-side 
shape for 4 < p% 19 < 6 GeV/c in the 0-12% most central 
collisions. Similar results were obtained for 3 < p^ 19 < 4 
GeV/c (not shown). 

The peak positions A in Fig. [5] show a slow increase 
with p°^ soc for the fits with the Symmetric Gaussian 
form. This functional form alone however, does not pro- 
vide a good description of the away-side shape for larger 
p^ 19 and p9j? soc . When an away-side contribution at 
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FIG. 5: (Color online) The angle A between the away-side 
peaks and A(p = n for 4 < p^,' a < 6 GeV/c in 0-12% cen- 
tral Au+Au collisions as a function of p^ saoc . Three different 
parametrisations of the away-side peak shape were used (see 
text). The lines show the systematic uncertainty from vi vari- 
ation while the errors are statistical errors from the fit. 



A<fi = 7r is included (2+1 Gauss and and 2 Gauss+Cosine 
shapes), the peak position A is close to 1.2 and approxi- 
mately independent of p T ssoc . This observation is quali- 
tatively consistent with predictions for a Mach Cone de- 
veloping in the hot and dense medium of the early stage 
of the collision [15| and with existing results on three- 
particle azimuthal correlations (39j . 



C. Associated particle spectra 

Figure [6] shows the integrated yield in the near-side 
peak (|A0| < 0.9) and away-side (|A</>| > 0.9) as a func- 
tion of p T ssoc in central Au+Au collisions and d+Au col- 
lisions for four ptf 19 intervals. The lower panels of the 
figures show the ratio of the associated yields in central 
Au+Au collisions and d+Au collisions. 

For d+Au collisions, the associated yield clearly in- 
creases with increasing p^f 19 and the slope decreases with 
ptrig ^ Thggg trends are expected from parton fragmen- 
tation, where the larger p^ lg selects larger underlying 
parton energies, thus increasing the multiplicity in the 
jet and leading to a harder fragmentation. 

The lower panels of Fig. [5] show that the ratio of the 
yields in Au+Au and d+Au is decreasing with p T ssoc , 
indicating that the fragmentation is softened due to in- 
medium energy loss. A softer fragmentation also implies 
that a trigger particle of given momentum selects differ- 
ent parton energies in Au+Au collisions than in d+Au 
collisions, which could explain some of the enhancement 
of associated yield at lower p T ssoc in Au+Au collisions. 
However, it should be noted that a large part of the in- 
creased yield at lower p T ssoc is at large Arj, associated 
with the ridge effect, and is therefore not necessarily from 



jet fragments [34(. It is also possible that other sources 
of particle production, such as parton coalescence and 
resonance decays, contribute at lower p^ 19 and may lead 
to different behaviour in d+Au and Au+Au. 



D. Azimuthal angle dependent mean-pr 

To further characterize the p^-dependence of associ- 
ated particle production, we perform an analysis of the 
inclusive mean-p^, (pt), of associated particles as a func- 
tion of Ac/>. The azimuthal distribution (pr)(Ac/>) is cal- 
culated by taking the ratio of the pr-sum distribution 
Pr (Ac/>) and the number distribution N(A<fi) 



(pr)(A$ =P T (Acj>)/N(A<f>). 



(1) 



The number distribution N(Acj)) = l/N trig dN/dA<j), as 
shown in Fig. [TJ while the pr-sum distribution Pr(A0) 
is formed using the same procedure, but addding the 
(scalar) transverse momenta as weights in the azimuthal 
distribution. 

To illustrate this method, Fig. [7] shows the inclusive 
distribution (p r )(A</>) for 0.25 < p T ssoc <4.0 GcV/c and 
three different ptf 19 selections for 0-12% central collisions. 
On the near-side, a clear increase of (pt) with p^ 19 is visi- 
ble, while the away-side {pr) distributions show a smaller 
dependence on p^ 19 . The lines in Fig. [7] show the back- 
ground. The elliptic flow of the background is calculated 



as a weighted average of v 2 



trig 



The difference be- 



tween the pT-weighted average (v2} PT , which is used to 
subtract the background in the pr-weighted distribution, 
and the number-weighted {v2)n, gives rise to the flow 
modulation of the background shown in the figure. 

To calculate the (px) of associated hadrons, the un- 
corrected background is subtracted from both the pr~ 
weighted and number-weighted distributions before tak- 
ing the ratio: 



(pt)(A0) 



P T (A0) - B PT (l + 2(v 2 ) PT cos(A0)) 
N(A<P)~ B N (l + 2(v 2 ) N cos(A<j))) ' 



(2) 



where N(Ac/)) and Pr(A0) arc the same number- 
weighted and sum-pT distributions used in Eq. Q] the 
average (1)2) are defined above and B PT and Bn are back- 
ground normalizations which are determined using the 
ZYAM method separately for the number and sum-px 
distributions. 

Figure [5] shows the resulting (px) of associated hadrons 
as a function of A<p in the away-side region for differ- 
ent ccntrality selections. In the peripheral bins a peaked 
structure in (j>t) is found, similar to the results in d+Au 
collisions (open circles). With increasing centrality, the 
{pr) around A(f> = n becomes lower. For the most central 
bin, the results show a minimum at A<f> = tt for the two 
softer trigger selections. For the highest trigger selection 
6 < p^ 19 < 10 GcV/c, a similar shape is seen, but there 
may be a slight enhancement at A0 = ir even in the most 
central collisions. 
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FIG. 6: (Color online) Near-side (|A<^| < 0.9, left panels) and away-side (|A0| > 0.9, right panels) associated yield per trigger 
particle for various p% lg selections as a function of px ssoc - Results are shown for 0-12% central Au+Au collisions and for d+Au 
collisions. The bottom panels show the ratios of the per-trigger associated yields in central Au+Au and d+Au collisions. The 
error bars on the points indicate the statistical uncertainty, including the statistical uncertainty on the subtracted background 
level. The gray bands indicate the uncertainty from the elliptic flow modulation of the background. 



We further study the difference between (pt) in the 
range \A<j> — tt| < ^ (referred to as "core" in the follow- 
ing) and at -| < \A(f> — tt| < f (referred to as "cone" 
in the following). Figure [9] shows (pt) in these two an- 
gular ranges as a function of the collision ccntrality for 
two different trigger px ranges. The (pt) decreases with 
ccntrality approaching the inclusive (pt), a feature al- 
ready reported in [To| for associated hadrons in the en- 
tire away-side region (| A</> — n\ < 2.14). This reduction of 
(px) is likely due to interactions with the medium. The 
fact that (pt) approaches the (pt) for inclusive particle 
production in events with a trigger hadron (solid lines in 
Fig. [!)]) may indicate that the associated particles at low 
Pt approach thcrmalisation with the medium. It is also 
clear that the (pr) of the core decreases more rapidly 
than that of the cone hadrons, which suggests that there 
is significant transport of associated hadrons away from 
A(j) = 7T due to jet-medium interactions. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, a comprehensive study of ccntrality and 
PT-dcpcndcncc of azimuthal di-hadron correlations in 
Au+Au events is presented. We observe several strik- 
ing modifications of the correlation structure in Au+Au 
compared to a d+Au reference. Associated yields on the 
near- and away-side are enhanced at lower pt- On the 
near-side, the increase in yield is partly located at large 
pseudo-rapidity difference A77 (see [3J] for a more de- 
tailed study) and the yields approach the measurement 
in d+Au collisions at the highest p^ 19 . On the away- 
side, the associated hadron distribution is significantly 
broadened; in fact, it is broad enough that it is impos- 
sible to unambiguously separate jet-like yields from the 
underlying event. At higher pt, p^? soc > 2 GeV/c and 
p^P 9 > 6 GeV/c, the away-side shape is narrow, like in 
d+Au events. A large enhancement of the away-side yield 
at low pt is found, while at higher px a suppression is 
seen with respect to d+Au collisions. 
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FIG. 7: (Color online) Inclusive, azimuthal (pr) of associated 
hadrons between 0.25 and 4.0 GeV/c in all events containing 
various classes of trigger particles. The lines show the V2- 
modulated background for the different trigger ranges, with 
the colors corresponding to the data points. 



These results are qualitatively consistent with a soften- 
ing of jet fragmentation by in-medium energy loss, lead- 
ing to an increase of the underlying parton energy se- 
lected by a trigger particle at given p 1 ^ 19 . Some of the 
changes in the correlation shapes could then be due to 
fragmentation of radiated gluons. 

The strong broadening of the away-side shapes, how- 
ever, does not seem to fit naturally in a description of 
particle production from medium modified jet fragmen- 
tation. Several alternative mechanisms have been pro- 
posed that could give rise to these structures. These can 
be divided into two categories: collective and radiative 
phenomena. 

Radiative treatments QJ, [H, EO, EH, S3] focus on the 
angular distribution of gluons radiated by the parton 
propagating through the medium. A large opening angle 
between the parent parton and radiated gluons is ex- 
pected when kinematic constraints are imposed. Two 
simplified calculations of this effect have been published 
in the literature. One of these calculations gives results 
that are qualitatively consistent with the data (l3j , while 
the other calculation [l2j shows a much smaller effect. 
Neither of the two calculations includes full integration 
over the initial state kinematics and the medium den- 
sity development. Another radiative scenario involves 
Cherenkov radiation of gluons [10, El], EH , which would 
give rise to conical distributions. 

Both for large-angle medium-induced gluon radiation 
and for gluon Cherenkov radiation, the expectation is 
that the a way -side shape becomes narrower with increas- 
ing p°^ soc [l3|,Ell ■ This trend is not observed in the corre- 
lation data: using a few different functional forms for the 
away-side distributions, we found that peak-separation A 
is approximately independent of p9j? soc . 

Alternatively, one could imagine that the passage of 



high-pT partons excites sound waves in the medium. It 
has been sug gested that this may lead to Mach shock 
waves [14T - [l7| . Qualitatively, the observed constant sepa- 
ration between the away-side peak and the constant coni- 
cal emission angle from three-particle correlations [39[ are 
consistent with this explanation. The transition from a 
broad away-side structure at low px to a narrow structure 
at higher pt would then signal the change from away-side 
structures dominated by bulk particle production from 
the medium to a situation where jet-fragments dominate. 

A recent 3D hydrodynamical calculation which in- 
cludes local density fluctuations in the initial state also 
shows a broad away-side structure that may be consis- 
tent with the experimental di-hadron correlation data 
[29j . In this model, there is no explicit introduction of 
hard partons or jets; the correlation arises purely from 
the medium. At the moment, it is not clear whether 
this model will also generate the conical emission signal 
seen in three-particle correlation data (39|. A study of 
three-particle correlations in this model is ongoing [44| . 

In general, a number of different mechanisms, includ- 
ing fragmentation, radiative energy loss, bulk response 
and hadron formation by coalescence of constituent 
quarks, may contribute to the observed di-hadron cor- 
relation structures. Quantitative modeling of the differ- 
ent processes, including the azimuthal correlation of the 
trigger and associated hadrons with the reaction plane, 
is needed to further disentangle the observed signals and 
the background. 

Experimentally, more insight in the underlying pro- 
duction processes will be gained from di-hadron mea- 
surements with identified particles and with respect to 
the reaction plane. In addition, 7-jet measurements, and 
measurements with reconstructed jets are being pursued, 
which provide better control over the initial state kine- 
matics. 
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FIG. 9: (Color online) Mean transverse momentum (pt)(A</>) 
of associated hadrons in the range 0.25 < p^ ssoc < 4.0 GeV/c 
in the core (\A(p - ir\ < f ) and cone (| < |A0 - vr | < §) 
azimuthal regions for 3 < p% lg < 4 GeV/c and 4 < p% 19 < 6 
GeV/c as a function of number of participants. The lines 
indicate the inclusive (pr) in the same pr range for events 
with a trigger hadron. 
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Appendix: Background shapes 



In this Appendix, we show the associated hadron az- 
imuthal distributions before subtracting the flow back- 
ground. 

Figure [TU] shows the distributions for associated parti- 
cles between 1.0 and 2.5 GeV/c and for different ranges 
of trigger particle px, ranging from 2.5 — 3.0 GeV/c (left 
column) to 6 — 10 GeV/c (right column) and different 
centralities (rows) . Note the large increase of the combi- 
natorial background with centrality. A large background 
modulation due to elliptic flow is expected. 

Figure [TT] shows the distributions of associated charged 
particles with various p T ssoc and p^ lg selections for 0- 
12% central Au+Au collisions. For p T ssoc > 1.0 GeV/c 
(lower three rows), the value of 1 / N tr i g dN / dA<f> at the 
minimum depends mostly on p^ soc , as expected for un- 
corrected background. In the upper row however, with 
0.5 < p T ssoc < 1.0 GeV/c, we observe a significant depen- 
dence of 1 / NtrigdN I dA<J) at the minimum on p^ 19 . The 
value at the minimum increases for lower p% %9 , which is 
qualitatively consistent with a centrality bias combined 
with the fact that the probability to find more than one 
trigger particle per event is sizable for the lower p^ 9 

selections and decreases with increasing p^ 19 . 

The background normalization (B) and the elliptic 
flow (v°; ssoc }(v2 n9 ) that were used to subtract the back- 
ground in Figs [2] and [3] are reported in Tables fl] and HU 
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FIG. 10: (Color online) Azimuthal angle difference distributions for associated particles with pr between 1.0 and 2.5 GeV/c 
and for different ranges of trigger particle pr, ranging from 2.5 — 3.0 GeV/c (left column) to 6 — 10 GeV/c (right column). 
Results are shown for Au+Au collisions with different centrality (rows). The line and the grey band show the elliptic flow 
modulated background that was subtracted to obtain Fig. [21 



15 




FIG. 11: (Color online) Azimuthal angle difference distributions for different 19 (columns) and p%? soc (rows) in 0-12% central 
Au+Au collisions. The line and the grey band show the elliptic flow modulated background that was subtracted to obtain Fig. 

El 
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PV 9 B|A,|<2.0 {V2 SSOC ){4" 3 ) Pt" 9 5|A„|<2.0 (l*""") «**> 

( GeV/c ) ( 10' 3 ) ( GeV/c ) ( IP' 3 ) 







60-80% 






0.5 <p^ ssoc <1.0 GeV/c 






2.5 


- 3.0 


1.220 ± 0.002 


20.6 ± 4.6 


2.5 


- 3.0 


59.833 ± 0.002 


2.9 


± 


0.8 


3.0 


- 4.0 


1.214 ± 0.003 


19.4 ± 4.6 


3.0 


- 4.0 


59.638 ± 0.003 


2.8 


± 


0.8 


4.0 


- 6.0 


1.197 ± 0.008 


15.7 ± 4.1 


4.0 


- 6.0 


59.366 ± 0.009 


2.6 


± 


0.8 


6.0 ■ 


- 10.0 


1.139 ± 0.026 


8.67 ± 2.7 


6.0 ■ 


- 10.0 


59.235 ± 0.034 


1.6 


± 


0.6 






40-60% 






1 


.0 <p^ aoc <1.5 GeV/c 






2.5 


- 3.0 


3.846 ± 0.002 


28.6 ± 4.4 


2.5 


- 3.0 


17.286 ± 0.001 


5.0 


± 


0.1 


3.0 


- 4.0 


3.835 ± 0.003 


27.5 ± 4.6 


3.0 


- 4.0 


17.238 ± 0.002 


4.9 


± 


0.1 


4.0 


- 6.0 


3.820 ± 0.008 


22.9 ± 4.5 


4.0 


- 6.0 


17.182 ± 0.005 


4.4 


± 


0.1 


6.0 ■ 


- 10.0 


3.76 ± 0.03 


13.0 ± 3.2 


6.0 ■ 


- 10.0 


17.154 ± 0.018 


2.7 


± 


0.1 






20-40% 






1 


.5 < p?r aaoc <2.5 GeV/c 






2.5 


- 3.0 


9.522 ± 0.002 


25.9 ± 3.0 


2.5 


- 3.0 


6.245 ± 0.001 


6.8 


± 


1.9 


3.0 


- 4.0 


9.494 ± 0.003 


25.4 ± 3.3 


3.0 


- 4.0 


6.230 ± 0.001 


6.6 


± 


2.0 


4.0 


- 6.0 


9.466 ± 0.008 


21.8 ± 3.4 


4.0 


- 6.0 


6.213 ± 0.003 


6.0 


± 


1.9 


6.0 ■ 


- 10.0 


9.515 ± 0.031 


12.8 ± 2.7 


6.0 ■ 


- 10.0 


6.211 ± 0.011 


3.7 


± 


1.4 






0-12% 






2 


.5 < p^ asoc <4.0 GeV/c 






2.5 


- 3.0 


23.531 ± 0.001 


5.5 ± 1.5 


2.5 


- 3.0 


0.5210 ± 0.0002 


(S.l 


± 


2.4 


3.0 


- 4.0 


23.469 ± 0.002 


5.4 ± 1.6 


3.0 


- 4.0 


0.5212 ± 0.0003 


8.0 


± 


2.5 


4.0 


- 6.0 


23.395 ± 0.005 


4.9 ± 1.5 


4.0 


- 6.0 


0.5207 ± 0.0008 


7.2 


± 


2.4 


6.0 - 


- 10.0 


23.365 ± 0.021 


3.0 ± 1.1 


6.0 


- 10.0 


0.5199 ± 0.0030 


4.5 


± 


1.7 



TABLE I: The v 2 and the normalization values used for the TABLE II: The v 2 values and the normalization used for the 
background subtraction in Fig. H background subtraction in Fig. 



